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INTRODUCTION 

The ST/MST cl ear- air Doppler radar, or wind profiler, is becoming an 
important tool in observational meteorology because of its capability to 
remotely observe dynamic parameters of the atmosphere. However, there are 
difficulties, which have long been recognized in work with precipitation 
sensitive Doppler radars, in transforming the observed radial velocities into 
meteorological wind components. As WALDTEUFEL and CORBIN (1979) put it ”... 
the problem lies in the fact that one would like to know as much as possible 
about the wind vector field, whereas a single-Doppler radar yields only one 
wind component, the radial. Every attempt to gain some knowledge of the vector 
field, therefore, must compound the data field with additional hypotheses or 
simplifications." In this paper, we review how this problem has been treated 
in the past, and recast some of the analysis to a form more suited to the high 
elevation angle, fixed beam ST radar profiling techniques. We then examine the 
diagnostic abilities of a number of fixed bean configurations with reference to 
a linear wind field. The results, in conjunction with other work which treats 
problems such as the effects of finite sample volumes in the presence of 
nonhomogene ous atmospheric reflectivity (e.g., KOSCIELNY et al., 1984), have 
implications important to the design of both individual MST/ST radars and 
MST/ST radar networks. 

BACKGROUND 

The use of pulse Doppler radar to measure horizontal winds seems to have 
been first suggested in the literature by PRGBERT- JONES (1960) in connection 
with a 3 -cm precipitation sensitive system. LHERMITTE and ATLAS (1961) 
suggested the VAD (velocity azimuth display) method of retrieving the mean 
horizontal velocity from radial velocity data taken in horizontal circles 
centered along the vertical of the radar site. GATON (1963) and BROWNING and 
WEXLER (1968) confined the analysis to a horizontal plane to deduce mean 
convergence, stretching and shearing deformation as well as the mean horizontal 
velocities from VAD observations. All of these authors treated only the 
stratiform situation, where 9w/3x = 3w/3y = 0.0. 

EASTERBROOK (1975) examined the processing of data in a conical sector, 
wherein values are estimated at a point not centered on the radar. In this 
case, although five parameters of the wind field can be extracted, the mean 
horizontal velocity is contaminated by vorticity and cannot be independently 
determined. 

The papers by WOODMAN and GUILLEN (1974) describing results using the 
Jicamarca incoherent- scatter radar, and GREEN et al. (1975) describing the 
results from the prototype VHF ST wind profiling radar at Sunset, Colorado, 
mark the arrival in the literature of VHF clear-air Doppler radars as important 
wind profile measurement tools. These radars are generally f ixed-multibeam 
systems, utilizing large elevation angles consistent with wind height- profiling 
through the troposphere and into the lower stratosphere, as opposed to the 
conventional meteorological radars, which usually have a single rotating beam, 
and utilize low zenith angles consistent with measuring winds at low heights. 
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Because precipitation sensitive Doppler radars work under conditions of 
relatively high signal to noise, a practical radar may utilize a rotatable dish 
antenna. Then it is natural to apply data-redundant least squares techniques 
taking advantage of the easily produced circular symmetry, such as the VAD 
technique (BROWNING and WEXLER. 1969). Although dear-air Doppler radars must 
work with the much .poorer signal- to-noise ratio provided by echoes from 
irregularities in the refractive index of the air itself, special high 
performance clear^air radars are able to create sufficient powers aperture by 
using very high power to allow use of a rotatable dish antenna. Under these 
conditions, VAD techniques can be used in clear air, as demonstrated by 
PETERSON and BALSLEY (1979). They used the Chatanika 23 cm Doppler radar to 
compare the accuracy of VAD, VED (velocity elevation display), and direct 
vertical probing in measuring w, the vertical wind component. They found 
direct measurement with the vertical beam superior to the other two methods. 

With typical ST radars, however, the power^aperture requirements couple 
with economic considerations to dictate large, immobile antennas for economical 
systems. It is not practical to steer these antennas by physically moving 
them, and though electronic beam steering methods are available (GREEN et al., 
1984, CLARK and GREEN, 1984. FUKAO et al.. 1985), they have not yet been 
generally applied. Thus, each beam position is expensive to implement, so that 
economical ST systems are designed to obtain the minimum amount of data to 
measure the wind components with as little bias as possible. 

CLARK et al. (1985) show the significant bias reduction obtained in u, v 
estimates over mountainous terrain obtained by adding a vertical beam to 
measure w. It will be apparent from the next section that the addition of 
additional beams can reduce the bias even further when stratiform conditions 
prevail. Hcwever, as discovered in the previous work with precipitation 
sensitive radars, it is not possible to eliminate bias completely with single- 
Doppler radar techniques. The following discussion will try to clarify the 
nature of this bias, and serve as an aid in design of economical ST radar 
systems. 

THE GENERAL SCALING EQUATION FOR A LINEAR WIND FIELD 

We adopt the usual Cartesian coordinate system x, y, z, representing 
distances to the east, north and zenith, respectively, with origin at the 
center of the radar antenna (Figure 1) . The primary assumption applied to the 
wind field is linearity in the region about (x , y , z ) , the point in 
space at which the flow parameters are to be determined. The vector distance 
to a measuring volume from the antenna is defined as 

r = xi + yj + zk 

where i, j, and k are the usual unit vectors. The wind vector at r is 
v = ui + vj + wk 
where it is assumed 

u = u o + (x-x o )u x + (y-y Q )u y + (z-z o )u z 

V = v o + (x-x o )v x + (y-y 0 )v y + (z-z Q )v z (1) 

W = w o + (x-x o )w x + (y-y 0 )w y + (z-z o )w z 

Th e x, y, z subscripts denote partial differentiation (e.g,, u = 3u/9x); 

U Q . v Q , and w q are the primary parameters we wish to measure, Eeing the 

values of u, v, and w at the point (x y z ). 

o ; o o 
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Figure 1. 


Often at this point in the analysis, a transformation to polar coordinates 
convenient to the operation of a radar is made. Instead, we will stay in 
Cartesian coordinates to facilitate the analysis. If we intend to make 
measurements of wind velocity at points (x q , y Q * z q ) for large ranges not 
centered on the radar antenna, it would be necessary to consider the curvature 
of the earth (DOVIAK and ZRNIC, 1984, p 261). Here we confine analysis to 
profile measurement directly above the antenna, where z is identical to height. 

The mean radial velocity v observed within a sample volume at a given r 
is related to the wind vector v ^here by the dot product 

rv = v*r = ux+vy+wz 
r 

Inserting the linear flow relations of equation (1) into equation (2), we 
get the general relation relating the parameters of the linear wind field, 
which we would like to measure, to the observed v^s: 

rv r = (u Q + u x (x-x o ) + u y (y-y o ) + u z (z-z q ))x 

+ (v Q + v x (x-x Q ) + v y (y-y Q ) + v z (z-z Q ))y 

+ (w o + w x (x - x o } + w y <y-yo> + w z^- z o ))z 


( 3 ) 
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Equation (3) is an arrangement of the general scaling equation 
convenient for developing sample location strategies. Much of the rest of the 
paper is concerned with describing some of these strategies. First, however, 
we will examine a second arrangement of the scaling equation in which general 
characteristics of solution are more apparent: 

rv r = *[u 0 -u x x 0 -u y y 0 -u 2 z 0 ] 

+ y tv o -v x x o -v y y o -v z z o^ 

+ ztw <f w x x o-Vo' W z z o ] 

+ * 2 Eu x ] <*> 

+ y 2 [v y ] 

+ z 2 [w ] 
z 

+ xy[u y + v x ] 

+ xz [u + w ] 

Z X 

+ yz[v + w ] 

J z y J 

This is a linear equation with nine parameters (the quantities in brackets), 
which may be solved from simul taneous observation of v at nine spatially 
independent locations. Of these nine parameters, only three correspond 
directly to the twelve parameters that specify the linear wind field. The 
other six are couples of these parameters, so that, in general, only linear 
combinations may be found. Specifically, solution of equation (4) for any 
given position (x q , y Q , z ) provides: _^l)the parameters necessary to 
describe the divergence of the flow V- v = [u ] + [v ] + [w ] ; 2) the 
additional three parameters [u + v ] , [u + X w ], aXd [v i w ] 
necessary to specify the deformation of tie flow; and 3) Z paraiieters containing 
u^, v Q , and w q coupled inextricably with spatial derivatives. 


PROFILING: MEASUREMENT OF THE WIND FIELD ALONG THE VERTICAL 

In volumes directly above the radar we have that x = y = 0.0, so 
equation (3) becomes ° ° 

rv r = x[u o + xu x + yu y + (z-z^uj 

+ y[v o + xv x + yv y + (z-z o )v z ] 

+ z[w q + xw x + yw y + (z-z Q )w 2 ] 

Once again, this equation may be rearranged to give 


(5) 


rv r = *[u 0 - z o u z ] 

+ y[v o - z o v z ] 

+ z[w - z w ] 
o o z 

+ x 2 [u x ] 

+ y 2 [v y ] 


( 6 ) 


C- (e 



389 


+ z 2 [w ] 
z 

+ xy[u y + v x ] 

+ xz[u z + w x ] 

+ yz[v z + w y ] 

where the bracketed terms represent the nine parameters that may be evaluated. 
The difference between this equation and equation (4) is the disappearance of 
all but the vertical shear terms from the first three parameters, leaving u , 
v q , and w coupled only to u z , v z , and w^, respectively. From the ° 

last two terms we see that when the stratiform condition is present (i.e., w 
and w may be neglected), u and v may be evaluated, thus allowing x 

unbiased solution for u anS v . Z 

o o 

VERTICAL BEAM MEASUREMENT OF w and w 

o z 

Measurements made with a vertically pointing beam are particularly simple 
to analyze, since they only contain information on two of the linear flow 
parameters: w q and w z * Inspection of the scaling equation with x = y = 0 
shows that the observed v is identically w , the vertical velocity 
component at height z . furthermore, the values of v observed at two 
closely spaced heights, z and z , yield w , so that r 


and 



(7) 


w 

z 


Thus a vertically direct ed Deam allows unbiased determination of w 
in the presence of linear flow. 


o 


and 


OBLIQUE BEAMS 


Assuming that w q and w z have been determined using a vertically 
pointing beam, we may treat these values as knowns and move them to the left 
side of the equation. It is then convenient to define the quantity as 


V = rv 
r r 


z [w - z w ] 
o o z 


2 

z w 


which contains only known quantities. Then we can write equation (6) as 


V 

r 


x[u o " 2 o u z ] 

+ y[v o - 2 o V z ] 
+ x 2 [u x ] 

+ y 2 [v y ] 

+ xy [u + V ] 

y x 


(8) 
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+ xz[u + w ] 

Z X 

+ yz[v z + w y ] 


This is the general overhead scaling equation when the vertical velocity and 
vertical velocity shear are known or negligible. It has seven parameters, 
requiring seven spatially independent observations for solution. 


ORTHOGONAL VERTICAL PLANES 


If measurement of the deformation of the flow is not required we can 
require all beans to be either in the xz or yz vertical plane. Then either 
y = 0 or x = 0 for all beans, so that the fifth term of equation (8) dis- 
appears. Under these conditions, equation (8) divides into two independent 
equations. For y = 0 


V 

— - [u + z u 1 + x[u ] + z[u + W ] 
x o o z x z x 

and for x = 0 


(9) 


V 

J- = t V o + W 1 + y[ V + Zlv z + W y 1 

In each plane, we now have three unknown parameters, requiring three spatially 
independent observations for solution. Thus, considering both planes and the 
two vertical beam observations used to evaluate V , eight independent 
observations are needed altogether when using orthogonal vertical planes. This 
does not mean we need eight different beams, since observations may be 
independent if they are separated in range along a beam. This technique will 
be considered in the next section. 


SINGLE BEAM ANALYSIS 

It would be very efficient if equation (9) could be solved by utilizing 
observations in three sample volumes along a single beam. However, such a 
scheme cannot lead to complete solution since z and x are not independent, 
being related by 

z = x ctn0 

where 0 is the zenith angle of the beam. For example, in the xz plane 
substituting for z using this relation yields 


V 


r 


x [u — z u ] 
o o z 
o 

+ x [u^ + ctn0 (u z + w^)] 


( 10 ) 


The number of parameters is thus reduced to two. While u q is still biased by 
z u , the stratiform approximation no longer allows solution for u , 
w£iih is now coupled to both u and u • Thus, neither u q nor V*v may be 
measured in unbiased form in tSe presence of u^. c 


CONSTANT HEIGHT ANALYSIS 


Starting again with equation (5), the vertical shear terms may be 
eliminated if we require that all observations be for z = (i.e., constant 

height analysis), giving 
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w r = xt u 0 + xu x + yu y ] 

+ y[v Q + xv x + yv y ] ( 11 ) 

+ z [w + XW + yw ] 
o o x y 

Gathering like terms 


rv 

r 


z 

o 


[w o ] 


+ x[u o + z o w x ] 

+ x 2 [u x ] 

+ xy[u y + v x ] 

+ y [v o + z o w y ] 

+ y 2 [v ] 
y 

This is a linear equation with six unknowns. The bias to u and v is now 
the z Q w x and z w terms, respectively, so that once again tfie strati- 
form approximation allows solution for u and v . 


( 12 ) 


USE OF THE VERTICAL BEAM 


It is easy to see that for a vertically pointing beam (i.e., x = y = o). 


v = w (13) 

r o 

Thus, w q is measured directly, with no bias. 

USE OF ORTHOGONAL VERTICAL PLANES 


Restricting beam positions such that either x = 0 or y = 0 splits equation 

, 

iwu luue^eiiueu l c^ua Lxuub • 


w r - z o w o = xtu o + z o w x ] + x2[u x J (1*> 

and 

w t " z o w o = y [v o + z o w y 1 + y 2[v y ] 

Thus, two beams in the xz plane plus two beams in the yz plane plus a vertical 
beam yield [w q ] . [u q + z q w x ]» [v Q + z Q w y ] , [u x ] , and [v y ] . 

DISCUSSION AND CONCLUSIONS 


The ability of ST radars to measure vertical velocity is unmatched at the 
present time by any other technique. It is straightforward and without bias, 
even if the linear wind field assumption is poor. Only under fiscal restraints 
of the most severe kind, or for very special purposes, should this beam be left 
out of ST radar systems, since this parameter is not only of general interest, 
but is necessary to remove bias from u and v measurements (CLARK et al., 1985). 

The key parameters to uncoupling terms in the scaling equations are w 
an< * V Whenever the stratiform condition, which states that these two X 
parameters are negligible, is satisfied, a five-bean ST radar may determine 
unbiased (with reference to a linear wind field) values of u, v, and w for 
sample volumes directly above the radar. Furthermore, the divergence and 
partial deformation of the flow may be determined. 
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Three-beam systems can determine w and w z . but are unable to obtain u 

and v wind components uncontaminated by vertical sheer terms, even when the 

stratiform condition is satisfied* 
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